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HIGHER ORDER ACCURATE PARTIAL IMPLICITIZATION: 

AN UNCONDITIONALLY STABLE FOURTH -ORDER- AC CURATE 
EXPLICIT NUMERICAL TECHNIQUE 

Randolph A. Graves, Jr. 

Langley Research Center 

SUMMARY 

An unconditionally stable fourth-order-accurate explicit numerical technique is 
derived, based on the method of partial implicitization. The Von Neumann stability 
analysis demonstrates the unconditional linear stability. The order of the truncation 
error is deduced from the Taylor series expansions of the linearized difference equations 
and is verified by numerical solutions to Burgers' equation. For comparison, results are 
also presented for Burgers' equation using a second-order-accurate partial -implicitization 
scheme, as well as the fourth-order scheme of Kreiss. 

INTRODUCTION 

A great deal of effort has been expended in developing numerical methods to solve 
fluid-dynamic problems, and much of the recent work is compiled in reference 1. A 
large part of that work centered on obtaining accurate solutions to the fiuid-dynamic equa- 
tions for transient processes. In many cases where only the steady-state solution is 
desired, the transient analyses have been applied since the introduced time dependency 
changes the boundary -value problem to an initial-value problem. The initial-value prob- 
lem lends itself readily to solution by explicit methods which have highly desirable char- 
acteristics for vector processing computers (see ref. 2). The steady-state solutions 
obtained by marching the transient problem asymptotically to steady state can still be 
costly in terms of computer resources since the maximum marching step is generally 
limited by stability considerations. Some effort has gone into speedily up the transient 
phase, increasii^ stability, so that the steady state is reached more rapidly; reference 3 
typifies this approach. Published in reference 4 was a partial-implicitization technique, 
an unconditionally stable second-order-accurate explicit scheme, which was shown in 
reference 2 to be well suited for use on vector processing computers. The unconditional 
stability of the partial-implicitization technique allows for more rapid calculation of the 
steady state. 


The purpose of the present paper is to present a modification to the procedure of 
reference 4 which produces a fourth-order-accurate method with only small changes in 
the second-order-accurate scheme. Maintaining the same form of the solution as ref- 
erence 2 insures the method will still have the desirable features for vector processing 
computers. 

SYMBOLS 

b intercept 

E average error (see eq. (20)) 

G amplification factor 

h spatial step size, Arj 

L total number of finite-difference points 

m order of accuracy 

R(. cell Reynolds number, VqW v 

R,S coefficients in scheme 2 finite-difference relations 

r,s coefficients in scheme 1 finite -difference relations 

t time 

t transformed time 

U velocity 

Uq velocity variable, U - U 

U steady-state wave speed, ^ 

V amplitude factor 

X coordinate direction 

a parameter in scheme 3 (see eq. (23)) 
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?7 transformed coordinate 

6 phase angle 

V viscosity 

Subscripts: 
e exact 

j finite -difference nodal point index 

Superscript: 

N time step index 

The primes indicate differentiation with respect to the transformed coordinate tj. 

MATHEMATICAL DEVELOPMENT 


Model Equation 

Burgers' equation has been used by a number of authors (for example, refs. 5 to 7) 
as a model equation with which to test numerical techniques. The equation is 


9t 9x 8x2 


( 1 ) 


In this form, equation (1) represents a diffuse shock wave through the application of the 
following boundary conditions: 


U(x,t) = 1 (x — -") 

U(x,t) = 0 (x — +“) 

Since only the steady-state solution is desired, the following wave-oriented transformation 
is applied to equation (1): 

77 = X - Ut 

A 

t = t 
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where U is the steady-state wave speed. Equation (1) now becomes 
8 t ° 9??^ 


where Uq = U - U. The boundary conditions under the transformations become 

(77 - -°o) 

iv - +°°) 



( 2 ) 


(3) 


The exact steady-state solution to equation (2) subject to boundary conditions (3) is 


U(? 7 ) = ^(1 - tanh ^ 
2 \ 4i^; 


(4) 


Solutions to equation (4) for v = 1/8, 1/16, and 1/24 are given in figure 1 where the 
decrease in wave thickness with decreasing viscosity is readily apparent. 


Partial Implicitization 


Scheme 1 (second-order accurate).- The derivation of the second-order-accurate 
partial-imphcitization technique as given in reference 4 will be repeated here in order to 
provide a base for the development of the fourth-order scheme. In reference 4, three- 
point equally spaced central finite-difference relations of second-order accuracy were 
used to express the first and second spatial derivatives of equation ( 2 ) at the points j - 1 , 
j, and j + 1. These finite-difference relations are: 


f^\ - Uj-2 

\9n /.I 2 At? 

' J -1 


/™\ ^.1 - »i-i 

\ 977 2 A 77 


I 

977^/. At]^ I 


(5) 


/£U\ _ ^j+2 ' 

\ 877 / 2 A 77 

j +1 


Aj\ ^1.2 - 2U.H1 ^ Uj 

8772/.^^ A772 
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Using the finite-difference relations (5) and expressing the time derivative as a simple 
backward difference, 

8U ^ 

8t At 

equation (2) can be written as 


. sluf.V 


+ (1 + 2s)u: 


N+l 

j 


- s - r 


]+l ] 


( 6 ) 


Uq -j AiL * V 

where r^ = — and s = Equation (6) is written at points j - 1, j, and j + 1; 

J 2 Ary Ary2 

however, the variables at points j - 2 and j + 2 are considered explicitly. (The 
explicit consideration of the variables at points j - 2 and j + 2 does produce an incon- 
sistent technique; that is, the resulting partial-implicitization technique cannot be used 
on true transient problems. This is of no consequence herein where the interest is in 
formulating a method to rapidly achieve a steady-state solution.) 

The system of three simultaneous equations is of the form 

Point j - 1 


(1 + 






+ 



+ s U 


N 

j-2 


Point 3 


-irj + s)U^_+^ + (1 + 2s)uj^+^ 



N+l 

j+1 



Point j + 1 



+ s U 


jN+1 


+ (1 + 2s)U 


N+l 

j+1 



+ 




(7a) 


(7b) 


(7c) 


In order to obtain the solution at point j, use Cramer's rule (see ref. 8) which gives an 
equation that will be used at all interior points in the finite -difference mesh except for the 
two points immediately adjacent to the boundaries. At these points, system (7) is solved 

to obtain equations for Not® that the equations for ^.nd 
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are used only for the two points adjacent to the boundary. The solution usii^ Cramer's 
rule is obtained from 


where 




(1 + 2s) 


-(r^ H-s 


(1 + 2s) 


rj + s 


A 

B 

C 


- '■j-i) 

(1 + 2s) 

•(■■i+i * =) 


0 

-(s - rj) 
(1 + 2s) 


f-i 

(1 + 2s) 


A = 




B = uf 


c = 


- ■'i.iK.a 


( 8 ) 


The result is 


= D 


(1 t 2s)uf + (s - rj)u^j + (s - 



+ 


(rj + s)u 


N 

j-1 





( 9 ) 


where 


D = 

(1 + 2s)^ 


1 

(s - ^ (s - + s) 
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Equation (9) now involves five points from the previous time level only; hence, equa- 
tion (9) is an explicit equation obtained from the partial implicitization of the difference 
form of the governing equation. In reference 4 this form was shown to be unconditionally 
stable for one dimension, and a similar form was shown to be unconditionally stable in 
the two-dimensional solutions of reference 2. In the present paper this partial- 
implicitization technique is referred to as scheme 1. 

Scheme 2 (fourth-o rder accurate).- In the following development, five-point finite- 
difference relations are used to express the first and second derivatives of equation (2) 
at the points j - 1, j, and j + 1. The use of five-point relations produce an overall 
accuracy of fourth order. These finite -difference relations are: 






-3Uj_2 - lOUj.i + 18Uj - euj^i . Uj ^2 


20 


f — 1 = 

1 r 

1,977 /. 

12 A77 


'j-2 ■ ■ ^j+2 


30 




-Uj.2 + Wj.i - 18Uj . lOUj^i . 3 Uj^2 


20 


9 2m 


V9j72/._^ 12 A772 


IIU ..2 - 20U..1 ^ - U.^2 


V 


M 1 

^977^/3 ^2 A772 L 


■U3.2 + 16U..1 - 30U3 + - u.^2 


-hluVi 

90 


'"92u\ 


1 


I }+l 12 At7^ L 


-Uj-2 + 4 U ..1 + 6 U 3 - 20U^^i llU.^2 


niuvi 

12 


( 10 ) 


Following a procedure similar to that used in the scheme 1 development, the finite- 
difference relations (10) with the time derivative expressed as a simple backward differ- 


ence 


9t 


U 


N+1 


U' 


N 


At 


-, equation (2) can be written at points j - 1, j, and j + 1. In 


scheme 1 the variables at points j - 2 and j + 2 were treated explicitly; similarly, 
the same procedure is followed here. However, it was found that the resulting technique 
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did not have unconditional stability; hence, some of the values of U at the points j - 1, 
j, and j + 1 had also to be treated explicitly in order to obtain unconditional stability. 
Again, as in scheme 1 the explicit consideration of the variables at the points j - 2, 
j - 1, j + 1, and j + 2 produces an inconsistent technique which is accurate only for the 
condition of steady state. 

U • At ^ 

When Ri = — and S = — ■ the system of three simultaneous equations 

) 12 A, ,2 4,2 

is of the form 
Point j - 1 

(1 . 20 S)u”V + (iSRj.j - 6 s)uN*‘ = uf.i e Rj_i(3uf_2 ' uf*2) 

* s(liuf .2 - uNj) r uN i(6Rj.i r 4s) 

+ (11a) 

Point j 

-^8Rj + i6sju?^Y * W + 30S)u|'+^ t (sRj - + *'i(''j+2 ■ ^ 1 ^ 2 ) 

- s(‘'jr2 * "i-2) 

Point j + 1 


-fl8R 


j+1 


+ 6S 




+ (1 + 20S)uf^/ = uNj ; Rj^l uf-2 - 3<2J ^ 


<1<2 


U 


N 

j-2 


<-ih.i - 


(11c) 


Use of Cramer's rule, as in the scheme 1 procedure, gives an equation which is used at 
all interior points in the finite-difference mesh, except for the two points immediately 
adjacent to the boundaries. At these points, system (11) is solved to obtain equations for 

and Note that the equations for and are used only for the 

two points adjacent to the boundary. The solution using Cramer's rule is obtained from 
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where 




(1 + 20S) 
■(8Rj + 16s) 
0 


A 

B 


(1 + 20S) 
■^8Rj + 16S) 


(iSRj.i - 6S) 
(1 + 30S) 

-(l8Rj+l + 6S 


(8Rj - 16S) 
(1 + 20S) 

0 

(8Rj - 16S) 
(1 + 20S) 


( 12 ) 


A = uf.i . 3uf_2 - UN 2 ) ^ sfliuf .2 - UN 2 , 


+ uf^l6R^.l + 4s) + 10R..iUf_i 


B = uf .RJuN2-uf-2 -sK2-uf-2 


^ = U^..R.Ju^-3uN,UsflluN.-uN_2 


]+l ]+ll ]-2 j+2j ]+2 ] 


- uf-lK.l - 4S) - lORj^iUf^i 


The result is 


uf ' . D /(I , 20S)|uf . - uf.j) - s(u"2 * ”" 2 ) ’ (™i ’ >®) 






uf-2)-<iKrt-«)-10R„iU! 


N 

j+1 


+ (8Rj + 16S 


uf.l - R,-i(3uf-2 - + s(uRf.j - u”j). u"i 


X (6R._^ + 4S^ + 


(13a) 
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where 


D = 


X 

(1 + 20S)(1 + SOS) - ^8Rj - 16s)(-18Rj+i - 6s) - (l8Rj_i - 6s)(-8Rj - 16s) 


(13b) 


Equation (13a) involves five points from the previous time level as in scheme 1; however, 
the difference relations in the present formulation involve all five of those points instead 
of only three. This technique is hereinafter referred to as scheme 2. 

Stability 

The Von Neumann stability analysis of equation (13a) can be performed by substi- 
tuting Fourier components of the form 

uj^ = V^e^^^ (14) 

The phase angle 0 is a function of the wave number kjj and the spacing h; that is, 


9 = kjj At] 

The phase angle varies in the range 


0 ^ e s 7T 


Substituting the appropriate forms of equation (14) into equation (13a) and assuming R 
and S to be constant results in 


yN+l _ 


-2S + 280s2 + 32r2) cos 20 + ^32S + 128s2 + 256 r 2] cos 9 


+ i(2R - 280RS) sin 20 - i(16R + 64RS) sin 0 + (1 + 20S) 


(15a) 


Looking first at equation (13b), 


D = 


1 

1 + SOS + 408s2 + 288r2 


(15b) 
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Since 0 < S < +<» and -«> < R < + 0 °, the denominator of equation (15b) is then always 

yN+1 

greater than zero and, thus, no singularity exists in equation (15a). Defining G = — ~ 

yN 

and noting that the Von Neumann condition for stability requires |g| = 1, equation (15a) 
becomes 


g| = D ^1 + 20S)2 + 2(1 + 20S)(-2S + 280s2 + 32 r2] cos 26 + 2(1 + 20S)(32S + 128s2 

+ 256r2) cos 6 + (-2S + 280s2 + 32 r2)^ cos2 20 + 2 (-2S + 280S2 + 32r 2)^32S + 128s2 


256r2) 


cos 6 cos 20 


\2 

+ ( 32 s + 128 s2 + 256r2) cos 2 0 + (2R - 280RS)^ sin2 20 

1/2 


+ (16R + 64RS)2 sin2 0 - 2^(2R - 280RS)(16R + 64RS) sin 20 sin 0 ]^ 


(16) 


Since equation (16) is rather complicated, consider first the simpler limits as in refer- 
ence 4. For S = 0 and R 0, equation (16) becomes 


G = 


1 + 4r2 + 256r2 + 64r2 cos 20 + 512R^ cos 0 + 1024R'* cos2 20 
+ 16384R^ cos 0 cos 20 + 65536R'* cos2 0 - 4 r2 cos2 20 

- 256r2 cos2 0 - 128r2 cos 0^1 - cos2 0j 


(17) 


1 + 288R'^ 


The maximum value of equation (17) occurs at 0 = 0, which results in 

. , (l + 576 r2 + 82944R^)^^^^ 

' ' 1 + 288r2 


Equation (16) satisfies the Von Neumann stability criterion in the limit S = 0 and R # 0. 
Similarly it can be shown that, when S # 0 and R = 0 (the maximum again occurs at 
0 = 0 ), 
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I I (l + loos + 33 les^ + 408008^ + 1664648^/^^ 

|G|=i '—‘1 

1 + 508 + 4088^ 

Equation (16) satisfies the Von Neumann stability criterion in both sets of limits, implyii^ 
equation (13a) is stable for all values of At in the limit conditions. It was determined 
by numerical evaluation that the maximum of equation (16) occurs at 6 = 0; thus, 


G 


^1 + 1008 + 33168^ + 408008^ + 1664648^ + 576r2 
+ 82944R'^ + 2880R^S + 235000r2s2) 

1 + 508 + 4088^ + 288R^ 


and, finally. 


|G|=1 

Equation (16) thus satisfies the stability criterion for all values of At and, hence, 
equation (13a) is an unconditionally stable solution for equation (2). 

Accuracy 

The formal accuracy of the linearized partial-implicitization technique, both 
schemes 1 and 2, is obtained from the steady-state form of equations (9) and (13a) by 
expanding each term of the partially implicit-difference scheme in a Taylor series 
expansion about Uj. The scheme 1 linearized steady-state counterpart to equation (9) is 

_ (s - r)^Uj^2 + (s + ^)\-2 
^ 2s 2 + 2r2 


8ubstituting the Taylor series expansions for Uj^2 ^j-2’ use of the linear- 

ized steady-state form of equation (2), and after algebraic simplification, the following 
result can be obtained: 


Uo 


£U 

dr] 




u:' 


V 

180 ^3 



(18) 
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The scheme 2 linearized steady-state counterpart to equation (13a) is 


(i40s2 + 140RS + 16r2)u._2 + (i28r 2 + 32RS + 64s2)u._j 

+ (i28r 2 - 32RS + 64s2)uj^^ + (l40s2 - 140RS + 16 r2)uj^2 

4083^ + 288R^ 


Following the previously outlined procedure for scheme 1 results in 


TT 9U 


8^U 


-34 . 4„ „V ^ 
U70 ^ 


o(h6) 


(19) 


Note that the third-order truncation errors from the second-derivative finite -difference 
relations cancel, thereby giving overall fourth-order accuracy for scheme 2. Thus, 
scheme 1 is of second-order accuracy whereas scheme 2 is fourth-order accurate. The 
increased accuracy of scheme 2 was accomplished with only a few changes to the scheme 1 
procedure. 

The formal accuracies given by equations (18) and (19) can be demonstrated numer- 
ically by using the following analysis. Defining the average error E as 


L-1 


E 




and assuming the average error to be a function of the spacing h = At] results in the 
following relation: 


L-1 


E 


= — - — )• U . - U. 

L - 2 ^ e,] j 


= bh 


m 


( 20 ) 


j=2 


where m is the order of the error term and b is considered constant. Taking the 
logarithm of both sides of equation (20) gives 


lnE = mlnh + lnb 


( 21 ) 


The value of m is determined by the slope of the plot of In E as a function of In h. 
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NUMERICAL RESULTS 


The solution to Burgers' equation was obtained over the interval -5 = t] = 5 with 
V = 1/8, 1/16, and 1/24 for a range of values of h. For the partial-implicitization tech- 

•A. ^ 

niques, the time step At = 1000 was used because this value exceeds the value of At 
necessary for convergence in the minimum number of iterations. For Kreiss' method 
(see ref. 9), At was taken to be a value which did not exceed the stability criterion. The 
solutions were considered to be converged when the following criteria were met: 




six 10-8 


The order of accuracy of schemes 1 and 2, as well as Kreiss' method included for com- 
parison purposes, was obtained from figures 2 to 4 where, from the analysis of equa- 
tion (21), the order m was taken as the slope of the corresponding curves. The partial- 
implicitization results indicate that the scheme 1 technique is of second-order accuracy 
and the scheme 2 technique is of fourth-order accuracy, thus verifying the formal accu- 
racies given by equations (18) and (19). The fourth-order Kreiss technique is more accu- 
rate in terms of the average error than the fourth-order scheme 2; but, the explicit nature 
of scheme 2 offers a number of advantages, particularly for vector processing computers. 
In tables 1 to 3 representative results are presented which reflect the relative accuracies 
seen in figures 2 to 4 for scheme 1, scheme 2, and Kreiss' method. Both of the fourth- 
order accurate techniques give considerably better results, compared with the exact solu- 
tion, than the second-order accurate method of scheme 1. This is particularly true for 
the smaller values of v and larger values of h which imply larger cell Reynolds 
numbers. 

From the solutions presented in tables 1 to 3, it can be observed that scheme 1 gen- 
erally produces a steeper wave than does the exact solution, whereas scheme 2 produces 
a more diffuse wave. This observation leads to the possibility of combining the difference 
relations of schemes 1 and 2 so that a wave more closely resembling the exact solution 
could be calculated. Assuming new finite-difference relations of the linearly combined 
form results in a solution referred to as scheme 3. Since 


Scheme 3 = o; (Scheme 1) + (1 - a) (Scheme 2) 

then the error terms in equations (18) and (19) involve a linearly. Next, setting a 
such that the derivative error terms through the fifth are zero results in the functional 
form 
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<y = 


( 22 ) 


^1 + 


where )3j and /32 are constants to be determined from numerical solutions. By vary- 
ing a over a range of values, the resulting average errors can be plotted as a function 
of a to determine that value of a which gives the minimum error for each Fig- 

ure 5 demonstrates just such a procedure. By calculating a number of a values cor- 
responding to the minimum errors, the constants in equation (22) can be evaluated. The 
best-fit results are given by - 7 and /32 = 0.88 so that 


a = 


R 


2 

c 


7 + 0.88R^ 


(23) 


Equation (23) is plotted in figure 6 along with values of a determined from error mini- 
mums. Equation (23) was used to determine the appropriate values of a and the test 
cases were rerun; these results are given in the last column of tables 1 to 3. It is 
immediately apparent that scheme 3 is always more accurate than schemes 1 and 2, 
particularly at larger cell Reynolds numbers. At these large cell Reynolds nupibers the 
combination of second-order-accurate and fourth-order-accurate schemes gives a fourth- 
order result, indicating that the lower order truncation errors may not be a true measure 
of what is occurring. The approximate fourth-order accuracy can be seen in figures 7 
to 9. 

The number of iterations to convergence varied from a minimum of 26 for scheme 1 
with i' = 1/24 and h = 0.1042 to a maximum of 1159 for scheme 3 with = 1/8 and 
h = 0.0521. For the conditions of tables 1 to 3, scheme 1 consistently took the fewest 
iterations whereas scheme 3 took the most. However, a more meaningful relationship is 
to compare the iterations to convergence for a given error level. From table 3, scheme 3 
took 114 iterations for h ^ 0.2083 to converge to an accuracy of 1.3860 x 10"^ and 
scheme 1 took 660 iterations for h = 0.0521 to converge to a comparable accuracy of 
1.2621 X 10““^. Scheme 2 took only 89 iterations for h = 0.2083 to converge to an accu- 
racy of 2.2395 X lO""^ which is in the range of the just-mentioned two cases. Thus, since 
schemes 2 and 3 do not require six times as many operations per step as does scheme 1, 
it appears that for a given level of accuracy, the higher order techniques are more 
efficient. 
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CONCLUDING REMARKS 


The second-order-accurate partial-implicitization numerical technique has been 
modified with little complication to achieve fourth-order accuracy yet retain the uncon- 
ditionally stable explicit feature of the method. The resulting fourth-order method still 
retains the desirable features for application to vector processing computers. In addition, 
an observation was made that at coarse grid spacings a linear combination of the second- 
and fourth-order schemes produces a more accurate result. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
July 8, 1975 
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TABLE 1.- SOLUTIONS TO BURGERS’ EQUATION FOR = 1/24 


(a) h = 0.2083; a = 0.500; Rc = 2.5 





U 



V 

Exact 

Scheme 1 

Scheme 2 

Kreiss’ 

method 

Scheme 3 

steady -state 

(second-order 

(fourth-order 

(linearly combined 


solution 

accurate) 

accurate) 

form) 

-1.2500 

1.00000 

1.00000 

1.00013 

0.99919 

1.00000 

-1.0417 

1.00000 

1.00000 

1.00056 

1.00162 

1.00008 

-.8333 

.99995 

.99985 

1.00153 

.99691 

1.00057 

-.6250 

.99945 

1.00137 

1.00090 

1.00590 

1.00208 

-.4167 

.99331 

.98780 

.98113 

.98180 

.99884 

-.2083 

.92414 

1.12500 

.86568 

.95055 

.92383 

0 

.50000 

.50000 

.50000 

.50000 

.50000 

E 


8.8688 X 10"3 

3.1649 X 10"3 

2.1574 x 10-3 ' 

3.9081 X 10“4 


(b) h = 0.1042; a = 0.186; R^ = 1.25 





u 



V 

Exact 

steady-state 

solution 

Scheme 1 
(second-order 
accurate) 

Scheme 2 
(fourth-order 
accurate) 

Kreiss' 

method 

Scheme 3 

(linearly combined 
form) 

-1.2500 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

-1.0417 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

-.8333 

.99995 

.99999 

.99998 

.99995 

.99998 

-.6250 

.99945 

.99985 

.99965 

.99944 

,99969 

-.4167 

.99331 

.99717 

.99402 

.99322 

.99473 

-.2083 

.92414 

.94944 

.91887 

.92415 

.92413 

0 

. 50000 

.50000 

.50000 

.50000 

.50000 

E 


1.6221 X 10‘3 

3.6103 X 10'^ 

5.5105 X 10-5 

1.7380 X 10"4 



(c) h = 0.0521; 

a = 0.0532; R, 

3 = 0.6250 





U 



V 

Exact 

steady-state 

solution 

Scheme 1 
(second-order 
accurate) 

Scheme 2 
(fourth-order 
accurate) 

Kreiss’ 

method 

Scheme 3 

(linearly combined 
form) 

-1.2500 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

-1.0417 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

-.8333 

.99995 

.99997 

.99996 

.99995 

.99996 

-.6250 

.99945 

.99957 

.99946 

.99945 

.99946 

-.4167 1 

.99331 

.99436 

.99335 

.99330 j 

.99340 

-.2083 

.92414 

.92999 

.92391 

.92415 

.92423 

0 1 

.50000 

.50000 , 

.50000 

.50000 

.50000 

E 


3.8364 X 10-4 

2.2888 X 10-5 

3.3731 X 10-5 

1.4873 X 10-5 




TABLE 2.' SOLUTIONS TO BURGERS' EQUATION FOR u = 1/16 


at 


+ Uo 


51 

dv 


a^u 


(a) h = 0.2083; oi = 0.294; Rc = 1.667 





u 



T? 

Exact 

steady -state 
solution 

Scheme 1 
(second-order 
accurate) 

Scheme 2 
(fourth -order 
accurate) 

Kreiss' 

method 

Scheme 3 

(linearly combined 
form) 

-1.2500 

0.99995 

1,00000 

1.00009 

0.99990 

1.00003 

-1.0417 

.99976 

1.00000 

1.00032 

.99987 

1.00018 

-.8333 

.99873 

.99993 

1.00018 

.99840 

1.00029 

-.6250 

.99331 

.99925 

.99442 

.99381 

.99701 

-.4167 

.96555 

.99180 

.95750 

.96359 

.96803 

-.2083 

.84113 

.91667 

.81639 

.84835 

.83507 

0 

.50000 

.50000 

.50000 

.50000 

.50000 

E 


4.6476 X 10-3 

1.5360 X 10-3 

4.3483 X 10-4 

6.3338 X 10-4 



(b) h = 0.1042: 

; a. = 0.0912; R 

C = 0.8333 





U 



n 

Exact 

steady -state 
solution 

Scheme 1 
(second-order 
accurate) 

Scheme 2 
(fourth -order 
accurate) 

Kreiss' 

method 

Scheme 3 

(linearly combined 
form) 

-1.2500 

0.99995 

0.99997 

0.99996 

0.99995 

0.99996 

-1.0417 

.99976 

.99986 

.99978 

.99976 

.99978 

-.8333 

j .99873 

.99917 

.99878 

.99872 

.99882 

-.6250 

.99331 

.99515 

.99344 

,99329 

.99360 

-.4167 

.96555 

.97206 

.96544 

.96553 

.96604 

-.2083 

.84113 

.85503 

,83925 

.84137 

.84059 

0 

.50000 

.50000 

.50000 

.50000 

.50000 

E 


1.0409 X 10-3 

1.1081 X 10-^ 

1.6336 X 10-3 

6.8578 X 10-3 


(c) h = 0.0521; o? = 0.0243; Rc = 0.4167 





U 



V 

Exact 

steady -state 
solution 

Scheme 1 
(second-order 
accurate) 

Scheme 2 
(fourth -order 
accurate) 

Kreiss' 

method 

Scheme 3 

(linearly combined 
form) 

-1.2500 

0.99995 

0.99996 

0.99995 

0.99995 

0.99995 

-1.0417 

.99976 

.99979 

.99976 

.99976 

.99976 

-.8333 

.99873 

.99885 

.99873 

.99873 

.99873 

-.6250 

.99331 

,99378 

.99331 

.99331 

.99333 

-.4167 

.96555 

.96716 

.96555 

.96555 

.96559 

-.2083 

.84113 

.84441 

.84102 

,84114 

.84110 

0 

.50000 

.50000 

.50000 

.50000 

.50000 

E 


2.5364 X 10“4 

1 

6.6485 X 10'^ 

9.8238 X 10"'^ 

4.4013 X 10-6 



TABLE 3.- SOLUTIONS TO BURGERS’ EQUATION FOR P = 1/8 




(a) h = 0.2083; 

a = 0.0912; R^ 

= 0.83333 





U 



V 

Exact 

steady -state 
solution 

Scheme 1 
(second-order 
accurate) 

Scheme 2 
(fourth -order 
accurate) 

Kreiss’ 

method 

Scheme 3 

(linearly combined 
form) 

-1.2500 

0.99331 

0.99515 

0.99344 

0.99329 

0.99360 

-1.0417 

.98473 

.98830 

.98484 

.98471 

.98517 

-.8333 

.96555 

.97206 

.96544 

.96553 

.96604 

-.5250 

.92414 

.93474 

.92333 

.92418 

.92434 

-.4167 

.84113 

.85503 

.83925 

.84137 

.84059 

-.2083 

.69706 

.70833 

.69506 

.69748 

.69613 

0 

.50000 

.50000 

.50000 

.50000 

.50000 

E 


2.1041 X 10“3 

2.2395 X 10-4 

3.3019 X 10“5 

1.3860 X 10^4 



(b) h = 0.1042 

; a = 0.0243; R 

c = 0.4167 





U 



V 

Exact 

steady-state ! 
solution 

Scheme 1 
(second-order 
accurate) 

Scheme 2 j 

(fourth -order 
accurate) 

Kreiss’ 

method 

Scheme 3 

(linearly combined 
form) 

-1.2500 

0.99331 

0.99378 

0.99331 

0.99331 

0.99333 

-1.0417 

.98473 

.98564 

.98474 

.98473 

.98476 

-.8333 

.96555 

,96716 

.96555 

.96555 

.96559 

-.6250 

.92414 

.92670 

.92411 

.92414 

.92416 

-.4167 

.84113 

.84441 

.84102 

.84114 

.84110 

-.2083 

.69706 

.69967 

.69694 

.69708 

.69700 

0 

.50000 

.50000 

.50000 

.50000 

.50000 

E 


5.0996 X 10“^ 

1.3367 X 10-5 

1.9743 X 10-5 

8.8487 X 10“5 


(c) h = 0.0521; a = 0.00617; Rc = 0.2083 


7 



U 



Exact 

steady -state 
solution 

Scheme 1 
(second-order 
accurate) 

Scheme 2 
(fourth-order 
accurate) 

Kreiss’ 

Scheme 3 

(linearly combined 
form) 

-1.2500 

0.99331 

0.99343 

0.99331 

0.99331 

0.99331 

-1.0417 

.98473 

.98496 

.98473 

.98473 

.98473 

-.8333 

.96555 

.96596 

.96555 

.96555 

.96556 

-.6250 

.92414 

.92478 

.92414 

.92414 

.92414 

-.4167 

.84113 

.84194 

.84112 

.84113 

.84113 

-.2083 ' 

.69706 

.69770 

.69705 

.69706 

.69706 

0 

.50000 

.50000 

.50000 

.50000 

.50000 

E 


1.2621 X 10-4 

8.1648 X lO***^ 

1.2053 X lO-*^ 

5.5739 X 10“*'^ 





Figure 2.- Average errors from solutions to Burgers’ equation for = 1/24. 
(The abbreviation PI indicates partial implicitization.) 
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Figure 3.- Average errors from solutions to Burgers’ equation for 
(The abbreviation PI indicates partial implicitization.) 


22 


= 1/16. 
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Figure 4.- Average errors from solutions to Burgers’ equation for v = 1/8. 
(The abbreviation PI indicates partial implicitization.) 
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Figure 7. 
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Average errors from solutions to Burgers’ equation for v = 1/24 
(The abbreviation PI indicates partial implicitization.) 




Figure 8.- Average errors from solutions to Burgers’ equation for p = 1/16. 
(The abbreviation PI indicates partial implicitization.) 
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Figure 9.- Average errors from solutions to Burgers' equation for v = 1/8. 
(The abbreviation PI indicates partial implicitization.) 



NASA-Langley, 1975 L- 10261 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON, D.C. 20546 


OFFICIAL BUSINESS 

PENALTY FOR PRIVATE USE S300 SPECIAL FOURTH-CLASS RATE 

BOOK 


POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 
451 



860 001 CL U G 730822 S00903DS 

DEPT OF THE AIR FCPCE 

AF ^vEAFCNS LABCRATCRY 

ATTN : TECHNICAL LIBRARY (SU L) 

KI RTLANC AFB NH 87 117 


POSTMASTER : 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 


"The aeronautical and space activities of the United States shall be 
conducted so as to contribute to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.^* 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 

Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. Also includes conference 
proceedings with either limited or unlimited 
distribution. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include final reports of major 
projects, monographs, data compilations, 
handbooks, sourcebooks, and special 
bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 
NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 




